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Abstract 

In this paper the geometric theory of separation of variables for 
time-independent Hamilton-Jacobi equation is extended to include the 
case of complex eigenvalues of a Killing tensor on pseudo-Riemannian 
manifolds. This task is performed without to complexify the manifold 
but just considering complex- valued functions on it. The simple for- 
malism introduced allows to extend in a very natural way the classical 
results on separation of variables (including Levi-Civita criterion and 
Stackel-Eisenhart theory) to the complex case. Orthogonal variables 
only are considered. 

1 Introduction 

The geometric theory of separation of variables for Hamilton-Jacobi 
equation on (pseudo)-Riemannian manifolds relates orthogonal sepa- 
rable coordinates with eigenvalues and eigenvectors of suitable Killing 
tensors. 

If Q is a real n-dimensional (pseudo)-Riemannian manifold with 
metric tensor g and real coordinates (q^), then any symmetric con- 
travariant two-tensor T defines on the cotangent bundle T*Q the func- 
tion T = \T^ u p^p v where p^ is the momentum conjugate to the coor- 
dinate q^. Two tensors Ti and T2 are said to be in involution or to 
commute in Schouten sense if the corresponding functions are in invo- 
lution with respect to the standard Poisson bracket on T*Q. A tensor 
K in involution with the metric is called a Killing tensor. Killing 
n-tensors play a crucial role in Riemmanian geometry because they 
encode the symmetries of the metric, they can equivalently be defined 
by the equation V ^K^ ... v \ = 0. 

In the study of separation theory it is usual to consider a single 
Killing two-tensor with pointwise distinct real eigenvalues and normal 
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independent eigenvectors (called a characteristic Killing tensor) or a 
family of n independent Killing tensors (including the metric) that 
commute algebraically and in Schouten sense (called a Killing-Stackel 
algebra). The main theorem 003113^ that geometrically characterize 
the separability of the Hamilton-Jacobi equation states: 

Theorem 1 The Hamilton-Jacobi equation associated to a natural 
geodesic Hamiltonian H = \g^ u PnPv is additively separable in the or- 
thogonal coordinates (q^) if and only if a characteristic Killing tensor 
K or, equivalently, a Killing-Stackel algebra exists, admitting the 1- 
forms dq^ as eigenforms. Being the coordinates orthogonal this is 
equivalent to ask that the vector fields are eigenvectors. The co- 
ordinate hypersurfaces coincide with the foliation orthogonal to these 
eigenvectors. 

The two approaches, based respectively on characteristic Killing ten- 
sors or on Killing- Sackel algebras, are interchangeable, and we will use 
one or the other according to which is more suitable to the situation. 

In pseudo-Riemannian manifolds (as the Minkowski space) real 
Killing tensors can have complex eigenvalues and eigenvectors. In 
previous papers the points of the manifold where this happens were 
considered as singular points, where no real separation coordinates 
were definable [Q. Separable complex coordinates were introduced 
in [2115] and subsequent papers on a complex Riemannian manifold. In 
this paper we introduce complex variables in a different way, avoiding 
the complexification of the original pseudo-Riemannian manifold. The 
complete integrals of the Hamilton-Jacobi equation found in this way 
turn out to be real-valued even if the variables used to determine 
them are complex- valued. This allow us to extend the definition of 
characteristic Killing tensor to pseudo-Riemannian manifolds without 
asking the eigenvalues to be real and do not require extension on 
complex manifold of the Hamilton-Jacobi theory. 

2 Complex variables from a Killing ten- 
sor 

Let K be a Killing tensor on the real re-dimensional pseudo-Riemannian 
manifold Q (with real coordinates (q^) and metric tensor g), such 
that n pointwise distinct functions and re functionally independent 
functions z^ are defined in a open subset of Q and satisfy (with no 
summation over fj,) the equations: 

Kdz>* = X^gdz^ . (1) 
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In other terms, K admits n independent exact eigenforms dz^ 1 , corre- 
sponding to the n pointwise distinct eigenvalues A^. We call this kind 
of Killing tensor a characteristic Killing tensor, dropping the request 
that the functions A^ (and as a consequence the functions z^) are 
real. Indeed, being the metric g not positively defined, the functions 
A^ and z^ can assume complex values. Because K and g are both 
symmetric with real coefficients, if z^ satisfies the relation (^Q) also 
its complex conjugate z^ satisfies the same relation. Hence, the func- 
tions z^ with a non-vanishing imaginary part are organized in pairs 
and the n complex- valued functions z^ define, through their real and 
imaginary parts, n real-valued coordinates. 

It is important to observe that also when K has complex eigenval- 
ues it is not necessary to complexify the manifold Q: the functions z^ 
are defined by using the real coordinates (q^), even if they have com- 
plex values. In this paper all the functions z^ are sometimes treated 
as independent, without taking in account the conjugation relation 
and using distinct indices for the function z^ and the function z^, also 
indicate by z^. 

This point of view is preferable to the use of the real and imaginary 
part of the functions z^ for at least two reason: because the simpler 
notation allow to handle the complex case in a way analogous to the 
ordinary real case and, more substantially, because the pairs of con- 
jugate variables can change from a point on Q to another. On a three 
dimensional manifold, for example, in each point at least one of the 
three function z^ is real-valued but any of them can assume complex 
values in some subset of Q. The use of different indices for a function 
and its complex conjugate allows to not specify which variables are 
real and which belong to a conjugate pair. In any case, in a neigh- 
borhood of any point of Q it is possible to split the functions z^ 1 in a 
subset of real functions and in a subset of conjugate pairs. The real 
ones will be denoted by latin indices: z a , while the conjugate pairs 
will be denoted by greek ones: z a and z a = z a . 

When the functions z^ assume complex values the definition of the 
derivation -J^ is no more obvious and some extra care is needed to 
use the z^ as a system of variables on the real manifold Q. Anyway, 
because of the conjugation relation, one must be careful to consider 
the z^ as coordinates. We will see that this abuse of terminology can 
lead to some incongruence, hence, we will call the functions z M simply 
"variables" . 



Considering the Jacobian matrix ( f^s" ) an d its inverse J one can 
define the, eventually complex- valued, vector fields 










(2) 



dz^ 
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Lemma 2 The vector fields satisfy the following properties: 

1. Z^q") = J"; 

2. Z^{z u ) = tf|; and if 0* / then Z^z**) = 0; 

3. ifz* = z^ then Z fl = Z fl ; 

4. [Z M ,Z„]=0; 



5. Z fl (W) = Z II (W). 

Proof. The first identity follows directly from the definition ©, ap- 
plying the vector field to the old coordinates (q u ). The second property 
holds because J is the inverse matrix of ^f^^ • The third is obtained 
from 



because = one has = and being Zp_ = J^-^p the thesis 
follows. To show that Z^ and Z v commute, a longer calculation is 
needed: from the definition of commutator one has 

/ d d \ d 

[Z ^ ZV] = \ J »df Jl ~ J »dq^ J ») W 

by differentiating with respect to q a the identity 

r^t - 8 X 



and setting in one case p = p, and in the other p = v, it follows that 

d 2 z x d 2 z x 
[Zp, Z u ](z x ) = Wq^p ~ J ° J Iq^ = • 

Finally, the last identity is proved observing that 



- dW dW 



Q q a H Q q a 



□ 



The vector field Z^ allows to differentiate a (possibly complex-valued) 
function / on Q with respect to the variable z^: 

l = 7m = r d]_ = d^dj_ 

For the variables z a and z a belonging to a conjugate pair the vector 
fields associated to their real and imaginary parts can be introduced. 
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Lemma 3 If z a = x a (q fi ) + iy a (q^) then the real-valued vector fields 

d d d 

+ 



dx a dz a dz® 

d . ( d d 



dy a \dz a dz' 

can be defined and they satisfies the relations -Q^x a = -^y a = 1 ; 

S*y a = w" x ° = °- Hence 

d 1 d id 



dz a 2 dx a 2 dy a ' 

d Id id 

+ 



(3) 



dz 51 2 dx a 2 dy° 



Proof. Being z a ^ z a the two vector fields Z a and Z& are different, 
thus both the vector fields and are well-defined. Applying 
these two vector fields to the functions x a = \{z a + z a ) and y a = 
^{z a — z a ) the desired relations are found. The expression of Z a and 
Za in terms of the vector fields and is a trivial consequence 
of the definition. □ 



The expressions (J2J) are sometimes called Wirtinger formal derivatives. 
Given z a and z@ both with non-vanishing imaginary part and (5 dif- 
ferent from a and a then, because x@ and y" depend only on z^ and 
z$, one has -^x? = ^/ = QE]. 

The formulas © are widely used in complex analysis. As instance 
if W = A(x fl ,y^) + iB(x^,y^) the equation -^W = is equivalent 
to Cauchy-Riemann equations. As an example of a misleading due 
to considering z^ and z^ as coordinates, let us consider the Hartogs 
Theorem ([H], Theorem 2.2.8): 

Theorem 4 Let f{z\ . . . , z n ) a complex function defined on a open 
subset of C, then f is holomorphic if and only if it is separately holo- 
morphic in each variables, i.e. if and only if 

df 

^ = 0, a = l,...,n. (4) 

It is obvious that W = Wi^) + W^z") is separately olomorphic, but 
if z u = z^ then does not follows that W is olomorphic. 

The efficiency of the introduced formalism becomes particularly 
evident in the definition of the momenta conjugate to the variables 
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through the formulas 

_dz^ ( 5 ) 

formally analogous to the canonical transformation associated to a 
change of coordinates in the manifold Q. If z^ is a complex variable 
also Pa is complex. Thanks to Lemma [2] one can define the conjugate 
momentum P v as: 

Jz h 

The canonical Poisson bracket on T*Q can be extended, by bilin- 
earity, to complex-valued functions and the transformation Q turns 
out to be canonical and to behave well under complex conjugation. 

The variables (P u , z^) are widely used in the literature on Hamil- 
tonian systems (e.g. see ^2])- The ma i n properties of these variables 
are recalled for sake of clearness. 



P» = J ^ = • (6) 



Lemma 5 The relations {P^, z u ] = 5%, {z^, z u } = and {P u , P v ] = 
hold. Therefore the previously defined variables {Pa, z^) are canoni- 
cal. 

Proof. The relations involving just the functions z^ are obvious, 
because they depends only on the (q^ 1 ). The other two relations are 

proved considering that the elements of the matrix ( &j ) depend only 



on the {q^) and using the standard argument: 

dq a dq a 
dq^dz^ _ 

- Q zl i Q q a - ft 
{Pp,P v } = i%iP^%Vr} 



dq a dq T dq a dq T 



dq T ,dq u ,dq u dq T ' 

dq a d dq T dq T d dq u 
dz^ Pr dq a dz u dz v PtJ dq T dz^ 
9<f 9 _ chf_ _d_ _ 

dz^ T dz»° dz" PrT dz» T ~ 



□ 
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Lemma 6 If = z M then Pp = P^. 

Proof. This property is immediately deduced taking the complex 

conjugate of the formula © and recalling that ||jr = |fjr and that 
the p^ are real- valued functions. □ 



The real and imaginary parts x a and y a of z a can be used to define 
the correspondent real- valued conjugate momenta. 

Lemma 7 If z a = x a + iy a has a non-vanishing imaginary part then, 
set z a = z, P a = P z , x a = x and y a = y, the complex momentum P z 
can be written as a function of the real momenta P x and P y , conjugate 
to the coordinates x and y, through the formula 

P Z = \(P X ~iPy) (7) 

Proof. If all the variables z^ with /j, different from a and a are left 
unchanged, the transformation 



x + iy z = x — iy 

dz dz 

dx z dx 
dz dz 

dy z dy 



Px — Pz ~\~ „ P z — Pz P Z 



p y = — p z + — p -z = KPz-Pz 



is canonical and defines the real conjugate momenta P x and P y . The 
formula (JTJ can be obtained, by solving the previous equations with 
respect to P z . □ 



In order to write in the variables (z^, P^) the Poisson bracket and the 
Hamilton equations, the vector fields 

d dz^ d 



dP^ dq» d Pu 

have to be introduced. As a consequence -jj- = J^-gp-- If z a has 
non-vanishing imaginary part then by setting z a = z = x + iy and 
being P a = P z = \{P X — iP y ) its conjugate momentum one can write 



d 

dP a 


d 
dP x 


d 


d 


d 


d 


dPa 


dP x 


l dPy 
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Lemma 8 The following properties hold: 
d d 



1. 



dPf, ' dP v 
d d 



dz^ ' dP u 



:0; 
0; 



3. if = then 



a 



a 



dP n op, 



Proof. The first relation follows because does not depend on p G , 
while the third is a consequence of Lemma |2j The second property is 



obtained by differentiating the identity 

dq T dz v 



S 1 ' 



z^ dq 7 

with respect to q G and observing that = o implies 

dq T d 2 z v 
dz^ dq T dq a 
Hence the thesis follows from 



0. 



a a 



dz» ' dP v 



dg T d 2 z v d 
dzt 1 dq T dq CT dp a 



□ 



Proposition 9 The Poisson bracket of two, possibly complex-valued, 
functions f(z^ 1 : P p ) and g(z^,P fl ) is given by 

df dg df dg 



{/,<?} 



dPp dz» 



dz» dP„ 



Hence, the Hamilton equations for an Hamiltonian H(z a ,P rT ) are 

dH 



Pu 



dP f 

_d 
dzv 



dH 



Proof. 



if, 9} 



df dg df dg 



dp a dq u dq a dp a 

dq a df dz v dg dz v df dq a 



dz^ dP^ dq a dz v dq u dz v dz^ dP^ 
df dg df dg 



dPp dz» 



dz^ dP„ 



□ 



Remark 10 Due to the conjugation relation between and z^ we 
have that the splitting of these complex Hamilton equations into real 
and imaginary parts yields two copies of the equivalent real Hamilton 
equations in real coordinates. This will happen again in the following 
and we denote the phenomenon as '"redundancy" '. 

The next Proposition is fundamental for separation of variables 
theory in complex variables: 

Proposition 11 The metric g and the characteristic Killing tensor 
K are simultaneously diagonalized in the variables (z' 1 , P„). In other 
words, the two associated quadratic functions 

H = \g' LV p^Pu , K = \K^ VilVv 
assumes the forms 



Proof. By the transformation formula (JSJ) in the new variables one 
has 

1 Bz u Bz T 1 Bz a Bz T 

2 y dq»dq v ' 2 dq>* Bq v 

but \ a g^B fJL z (T B v z T = K^d fl z a d v z T = X T g^B fl z a B u z T and thus (be- 
ing X a pointwise distinct functions) if a ^ r it follows g^ v B^z° B v z T = 
and K iiV d j j L z a d v z T = \ a g fJ,u B fl z a B u z T = 0. The thesis is obtained by 
setting, for the non-vanishing coefficients, 

Bz a Bz a Bz a Bz a 

g» v =g° K^ U ————=K a . 

BqV Bq v Bq^ Bq v 

□ 



Redundancy manifest itself here by the straightforward relations 
g ~« = g ° and K a = K°. 

Remark 12 Given a characteristic Killing tensor K^ , or the associ- 
ated Killing- Stackel algebra, the functions z^ are given by integrating 
either the eigenforms or the eigenvectors of K^ v or (as we will see) 
by the eigenvalues of a basis of the associated Killing- Stackel algebra 
through the fundamental functions f% c given in generalized to the 
complex case. 



9 



3 The Hamilton— Jacobi equation and 
its solutions 



The Hamilton-Jacobi equation for Hamiltonian H, written in complex 
variables, is obtained by setting 

p - dw (Z) 
^"^7- (8) 

in the expression of H. In this way a PDE with complex coefficients 
for a function W is found and in general its solutions will be complex- 
valued. The consistence of the substitution (jHJ) implies that Pp = P ^ 
if zp = Zn,- This fact is achieved only by using a real- valued solution 
of the Hamilton-Jacobi equation. 

Proposition 13 Let H(P^ l ,z fl ) be a time-independent Hamiltonian 
written in the canonical complex variables (P^,^). A solution W of 
the Hamilton- Jacobi equation associated to H : 

dW 

H(—,zn = h (9) 



has a constant imaginary part (and then an equivalent real solution 
exists ) if an only if is complex conjugate to g=£ . 

Proof. From Lemma 12 it follows that (set z^ = z^) 



dW _ Ju dW 

dW _ j u dW 
dzP- ~ ^Ihf 

Hence, being the matrix J invertible, the first equation is the complex 
conjugate of the second one if and only if 

dW dW 



dq u dq u 

that is if and only if the imaginary part of W is constant (and is 
equivalent to zero). □ 



The relation between the Hamilton-Jacobi equation with complex 
variables and the ordinary one can be clarified by using (when z a = 
z = x + iy has a non-vanishing imaginary part) the non orthogonal 
coordinates given by their real and imaginary parts. Being 

dVV_ldW_ i dW 
dz a ~ 2~dx~ ~ 2 dy ' 
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by using the formula (JJJ) one obtain the substitution employed in the 
ordinary Hamilton-Jacobi equation: 

_ dW _ dW 

dx ' y dy 

Definition 14 A complete integral of the equation $3) is a solution 
W(z fl ,c u ) of where c v are n real-valued parameters, satisfying 

A complete integral W (z 11 , c u ) is (additively) separated in the variables 
(z») if 

W{z», c) = W 1 {z\c u ) + --- + W n (z n , c). 
The definition © of implies 

\dz^dc v J 1 ' \dq a dc l 

then, because J is invertible, when is real the definition of complete 
integral in the complex case is equivalent to the usual one. It is im- 
portant to notice that by rewriting the function W in the coordinates 
(q^) one obtains different results if the variables {z^) are complex or 
real. The introduction of the complex variables allows to find a com- 
plete real integral of the Hamilton-Jacobi equation and then to apply 
the classical Jacobi method even where the eigenvalues of K are not 
real. 



det h^T = detJdet b^H (io) 



4 Separation of variables for Hamilton- 
Jacobi equation 

The Levi-Civita criterion for the complete additive separation of a 
Hamilton-Jacobi equation can be rewritten in our complex variables. 

Following Levi-Civita the Hamiltonians considered will be always 

dH 

non-degenerated, that is satisfying ~Qp~ 7^ for all the indices, but not 

necessarily natural. Moreover, it is also possible to consider a complex 
valued Hamiltonian. 

Theorem 15 The Hamilton- Jacobi equation for Hamiltonian H ad- 
mits a complete integral separated in the variables (z^), if and only if 
the Hamiltonian H satisfies the Levi-Civita equation 

d^Hd u Hd^ v H + dpHdyHd^H - d»Hd u Hd»H - d^Hd v Hd^H = 
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where fi ^ v, no summation over the repeated indices is understood 
and 

d dq v d ^ d dz^ d 
dz^ dq v ' Pa dq v dp u 



Proof. The proof is adapted from the Levi-Civita's one The 
Hamilton-Jacobi equation for the Hamiltonian H(P^, z^) is H{w 



h where P^ = = w fl (z u ). By differentiating this equation with re- 
spect to z v one obtains 



dH dH dw L 



+ 



dz v dP^ dz v 



0. 



The additive separation of variables is obtained if and only if de- 
pends only from z^. If Mf- ^ the complete separability is equivalent 
to the existence of a solution for each of the n PDE system 



Q Z » dH_ 
dP u 



dz u 



{y + m) 



(11) 



Each PDE system can be seen as a system for a function de- 
pending on n independent complex variables, but the variables z^ 
obtained from Q are not independent due to complex conjugation 
relations. Let us split the functions z^ in real ones and conjugate 
pairs: (z^) = (z a , z a , z a ). The system (jllj) can be rewritten using the 
Wirtinger formal derivatives (j3J) and, depending on the value of /j,, one 
of the following systems is obtained 



dw a _ p 
dz a — ri a 

8Wg 



dx a 

8Wg 



( dWg _ T> 

dWg 

dy a 

dxfi 
dw a 



dw a 



iR 







Rn 



8Wq 

dx a 

8Wg _ 

dy a ~ 

dWg _ Q 

dwg _ n 

dz a ~ U 



iRn 



(12) 







These systems admit a solution if and only if the following integrability 
conditions (with (3 different from a and a) respectively hold: 





dR a 


dR a 


_ dR a 




dz b 


dx a 


dy a 


dR a 


dR a 


dR a 


dR a 


dz a 


dxP 


" dyP 


dz* 


dR a 


dR a 


dR a 


dR Q 


dz a 


dxP 




dx a 



o. 

0. 
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Due to the linearity of the total derivative with respect to z^, the 
necessary and sufficient conditions for the existence of a solution of 
systems (fTTj) become: 

dR u dR u n dR u 

az u oz u ow v 

where no summation over the repeated indices is understood. Sub- 
stituting the expression of R^ from (jllj) the Levi-Civita condition is 
obtained. □ 



Proposition 16 If the Hamiltonian H is real and satisfies the Levi- 
Civita equations then the complete integral can be chosen to be real. 

Proof. The definition (jllj) of R^ proves that if H is real then R a are 
real and R a = R a . Hence, the first of the systems (j!2j) implies that 
w a can always be chosen real. The case of w a , instead, is more subtle: 
first one must observe that in the separated integral W the function 
W a always appears added to the function W a , then one must consider 
the function w a — w a that has the same imaginary part of w a + w a . 
The last two of the systems (j!2j) imply that 

d 



dx a 

d 



[W a - Wa) = R a - R a =0; 

(w a — Wa) = iRa — iRa = . 



dy° 

Hence the imaginary part of w a +w a = W a +W a is a constant that can 
be chosen to be zero and this means that also the imaginary part of 
W a + W a is a constant and it can be chosen to be zero. It follows that 
always exist a real complete integral of the Hamilton-Jacobi equation. 
□ 

Remark 17 The previous theorem holds for any system of canoni- 
cal variables (z^,Pu) independently from their association with eigen- 
forms of a Killing tensor and from the signature of the metric g, pro- 
vided the non real variables are pairwise conjugate. 

Theorem 18 The functions K {a) = \ ^ Af^P^ 2 eH = \Y, v 9 up ^ 
are in involution if and only if the Eisenhart conditions hold: 
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Proof. Using the Poisson bracket properties one obtains 
{K (a) ,H} = ~EAK(*),9"}P» 2 + Ev{K(*),P»}Pv9" 

= \E^U9^9 u }P, 2 Pu 2 + \Z»AP^9 U }^ P ^ 2 + 

jE^/.^TO 2 + \e,Ap, 2 ,Pu}xU^9 u p^ 

1 r r)\ I/ 1 

— 2^M" M M^Szm y dzv \9 r ii r v 
and the thesis follows from the identity of polynomials. □ 



The existence of a KS-algebra implies the equations (|!Hll5jl . The 
converse rises some monodromy problem in the definition of the func- 
tions \f x that will not be examined here. In the following, the ex- 
istence of a KS-algebra or of a characteristic Killing tensor will be 
understood, together with the associated complex variables. 

By writing the Hamilton-Jacobi equations for all the Hamiltonians 
Km = H, . . . , Kr n ) associated to a KS-algebra the following system 



of PDEs is obtained 



/ 9 1 
l 

(2) 



\1 „1 
A (2\9 



A (2)9 



\ ( (W 2 \ 

(d 2 W) 2 



A (n)9 I 



\ {d n Wf J 



( a \ 

\c n J 



(16) 



Denoting by S the inverse of the matrix in the left side of the equation, 
one has 

Y, v Sl9 v = 5l o = \ 



(17) 



Theorem 19 The Eisenhart equations 115)) hold if and only if the 
entries of the invertible matrix S satisfy ■JLtSJ^ = for fj, ^ v. 

Proof. The differentiation of the two equations (|17| ) gives, respec- 
tively 



1 



9 U ^F^Z + Ylu si gf? 9 V — o 
The second expression can be written as 



S;.^^(^ ct )^)-Z ; /^^ ct) ^^+^5^A^ ) ^^+E^A|' ct) ^^5J - 
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and subtracting the first expression (multiplied by A^ 1 J one has, with 



Therefore, if d^SJ, = for fi ^ v the second summation vanishes and 
the invertibility of S imply the Eisenhart equations. 

On the converse, taking the Eisenhart conditions as hypotheses for 
any fixed value of and r one has the system of n — 1 equations 



J2 [m m 



U 



in the n — 1 unknown ui^ T ^ = d^SJ, with v ^ fi. The determinant of 
the matrix 



M (m) 

la 1 



/ AL-A^ • • • A^-AiL Af+'-A^ • • • A?„,-A^ \ 



5 n / 



\2) "(2) 



(2) /v (2) *(2) ^(2) ^(2) ^(2) 



V A W" A (n) • • • \n) ~\n) A („) ~\n) ' ' ' A (n)" A (n) / 



IS 



X 1 \^ _1 _X<" \L l+1 —\t 1 \n _X^ 

A {2) A (2) ' ' ' A (2) A (2) A (2) A (2) ' ' ' A (2) A {2) 
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X 1 \^ \^ +1 —\> l X™ 

A {2) A (2) ' ' ' A (2) A (2) A (2) A (2) A (2) ' ' ' A (2) A (2) 



A/ \ — A/ \ * * * Aj \ — A { \ A/ \ — A { \ ' ' ' A/\ — At \ 
(n) (n) (n) (n) (n) (n) (n) (n) 



\1 i/i \/i+l \M \M \M \n \M 

A/ \ — Aj \ * * * A/ \ — A/ \ A/ \ A/ \ — A/ \ * * * A/\ — A/ \ 
(n) (n) (n) (n) (n) (n) (n) (n) (n) 



n< 



A {2) A {2) A (2) 



X 1 . . . X^ 1 . . . X n 

A («) A W A W 



ir 
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and it vanishes if and only if the determinant of S vanishes. The 
unique solution of the system (|16|) is then u^ r ^ = for all fj,, r, v. □ 

The matrix S, where the ^-th row depends only from is the 
complex generalization of the classical Stackel matrix. Taking the 
inverse of (|TH|) the separated ODE's 

8W dW u 

dz" dz" vCt: 

are found. Because the function W is real it is possible to apply the 
classical Jacobi method to the function W(z fM (q CT ),c u ) and solve the 
Hamilton equations. The complex variables are introduced only in 
order to find a separated complete integral of the Hamilton-Jacobi 
equation. Since this integral is found, being it a real function, every- 
thing can be done in real coordinates and a complex version of the 
Jacobi theorem is not needed. 



5 Examples 



5.1 Two-dimensional case 

Because complex variables with non-vanishing imaginary part always 
appear as conjugate pairs, the two dimensional case is particularly 
relevant and is a model for each pair of truly complex variables. 

Let us consider the particularly simple case of a two-dimensional 
manifold Q on which the function z 1 = z has everywhere non- vanishing 
imaginary part. Then the relation z 2 = z necessarily holds and the 
coordinates (x, y) are defined by 



The matrix J, in this case is 



x + iy 
x — iy 

dz 



2J -4i 




where 



dx dy 
dq 1 dq 2 



Jtf) = 

The vector field J| is given by 



dx dy 
dq 2 dq 1 



dz 
dq 1 
dz 

dq 2 



dz 
dq 1 

dz 

w 



d_ 

dz 



i 

2J 



dz d dz d 



dq 2 dq 1 dq 1 dq 2 
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and the vector fields e ^ can be written as 

d 1 / dy d dy d \ 

dx J \dq 2 dq 1 dq 1 dq 2 J 

d 1 / dx d dx d \ 

dy J \ dq 1 dq 2 dq 2 dq 1 J 



The link between old and new momenta is given by 

dz dz - 
P » = dq^ F + dq^ F 

that imply 

p -h < 18 > 

The Poisson bracket of P with z and z can be calculated directly 
from the formula (|18|) and the definition of J: 

dP dz dP dz 



{P,z} 



dp^ dq^ dq^ dp^ 
i ( dz dz dz dz 



{P, z} 



2 J \dq 2 dq 1 dq 1 dq 2 
i 2J _ 
2jT ~ 

dP dz dP dz 



dp ^ dqt 1 dqt 1 dp^ 
i ( dz dz dz dz 



2 J \dq 2 dq 1 dq 1 dq 7 
= 

as expected for a canonical transformation. 

Finally the Hamiltonians H and K associated to the metric and 
the Killing tensor are, in the x, y coordinates: 

~ 2 9 \dq^ x+ dq^ y )\dq" Px+ dq" Py 



" 2 9 dq»dq" Px+9 dq»dq" PxPv+ 2 9 dq» dq" Py 

but because x = l(z + z) e y = — z) one has 

v dx dx „j, / dz dz \ / dz dz \ 

9 dq^dq 1 ^ = 9 \dq^ + ^J \dq^ + dq l; ) =9 + 9 

^dq^dlf = ^ {dq^ + dq^) =%{9 ~ 9) 

, v dy_dy_ = _ f dz_ _ dz_\ ( dz_ _ dz_\ = 
9 dqVdq" 9 \dq» dq^ J \dq" dq" J " 9 9 
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that is 
H = 



1 / , 2 * ' - 1 ' s o 1 



(9 + g)pi + 4(5 - - I (5 + a)p'y = ^(gP 2 + gP 2 



1 - i 

K = ~{Xg + Xg)p 2 x + -{Xg - Xg)p x p y 



^{Xg + Xg)p 2 y = - 



XgP 2 + XgP 2 



Hence, the coordinates x and y (differently from the variables z and 
z) do not diagonalize neither the metric nor the Killing tensor. 

A complete separated integral of the Hamilton-Jacobi equation 
associated to H is a function W(z, z; h, k) such that, if c\ = h and 
c 2 = k, 

W{z,z;h,k) = W 1 (z;h,k) + W 2 (z;h,k) 
dW - dW" 



H P 



det 



dz 

f d 2 W 

\dz^dc u 



,Pz 



dz 



h 



The fact that K is a Killing tensor for the metric g is equivalent to 
the involutivity between the functions H and K, hence the equivalent 
Eisenhart conditions take the form 



dX u 
dz^ 



d 



(19) 



where no summation over repeated indices is assumed. The equations 
(|19J) correspond to the two (complex conjugate) systems 



f OX 

dz 







dX - dg 
9 dz- = {X - X) dz- 



dX 

dz 







-dX -dg 

dz dz 



or, equivalently, to the two systems 

d_ 

dz 
d 



X = 



k dz 



g(X - A) 







a 



dz 
d_[ 

dz 



X = 



g(X - A) 



(20) 







The two Hamilton-Jacobi equations for the two Hamiltonians H and 
K, in the variables z and z, give the PDE system 



9 



dW 

dz 



+ g 



dW 



dz 



h 



dW 

dz 



+ Xg 



dW 

dz 
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that is by solving with respect to the derivatives 




( m 2 ) - 

The systems (j2U|) imply that the first row of the matrix depends only 
from z, while the second one only from z, hence if and W2{z) 

satisfies 



8WA 2 


hX 


- k 


dz J 


<?(A 


-A) 


dW 2 \ 2 


hX 


-k 


dz J 


9(A 


-A) 



the function W(z, z) = Wi(z)+W2(z) is a real solution of the Hamilton- 
Jacobi equation for both the Hamiltonians H and K. 

It is important to observe that, being the complex conjugate 
of and being h and k real-valued, the function W2 satisfies the 
same equation of VFi(z) and so the solution of the Hamilton-Jacobi 
equation is found integrating just one equation. Again appears the 
phenomenon of redundancy. 

Moreover, if A is not constant from 1)2(^1 follows that A is a function 
of z only. By taking A = z it is possible to determinate a standard 
form for g: 

z — z 

where A is an arbitrary function of z and, consequently, olomorphic. 



5.2 First example: everywhere complex vari- 
ables 

Let consider, in the Minkowski plane with coordinates q 1 = x and 
q 2 = t, the metric and the Killing tensor 

co = (0 -1) («n = (! I) 

The corresponding Hamiltonian functions are 

H = \{pl-V 2 t) K = p xPt 

The equation K dz^ = A^g dz^ 1 relates to the eigenvalues dbi the func- 
tions z 1 = z = x + it and z 2 = z = x — it. The eigenvectors of K 
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have non- vanishing imaginary part everywhere the change of variables 
is invertible, thus no real separable coordinates can be found. 
The Jacobian matrix and its inverse J are 



fdz v 



then, the momenta p x and pt are given by 

p x = P + P p t = i(P-P). 
In the new variables the Hamiltonians become 

H = P 2 + P 2 K = i(P 2 - P 2 ) . 
The two Hamilton-Jacobi equations are 



1 1 

-i i 



dW 

dz 
dW 

dz 



+ 



dW 

dz 
dW 
dz 



h 



-ik 



and looking for a separated solution W = W\{z) + Wz(z) one has that 
W\ and W2 (taking h and k real) both satisfy the equation 



dl 



h — ik 



Set Q = h + ik (not belonging to the negative real semi-axis, to en- 
sure the monodromy of the square root) the solution of the previous 
equation is 



S = z 



that is (see §21.2 in 0): 



V2Re (zy/( 
V2(xRe^/C-tlm^ 



h+y / h 2 + k 2 + 



kt 



\fh + Vh 2 + k 2 

To find the conjugate momenta for h and k one can apply the 
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Jacobi theorem, being W, h and k real: 



9W 1 -Re(Wt 



dh |c|V2 

= K^( xRe ^ + tIm ^) 

Ok = ^77- = , ~ ,- Im (zy/C 



dk \£\y/2 

where 9j~ is constant while 9^ = r — tq. 

By solving with respect to x and t one has 

t = y/2 (OhBey/C - e k 1my/c) 
and the integral curves are (as expected) straight lines. 



5.3 Second example: parabolic coordinates 

Let us consider, in the Minkowski plane with coordinates q l = x and 
q 2 = t, the metric and the Killing tensor 

on =(j _°o <-) - ( ) 

The corresponding Hamiltonian functions are 

H = - (p 2 x - p 2 t ) K = -{p\ + p 2 t ) + ( Px + Pt ) (x Pt + t Px ) . 

The functions z\ , Z2 are found by solving the equation K dz^ = A^g dz M , 
that relates to the eigenvalues A M = t — x — (—1)^-^/1 + 2(x + 1) the 
functions = x — t— (—1)^-^/1 + 2(x + i). It is worth to observe that 
A 1 depends only from z 2 e vice versa. Moreover the two eigenvectors 
are real and distinct or complex conjugate depending on the sign of 
2(x + 1) + 1; the line x + t = — ^ consists of singular points where the 
eigenvectors coincide. Then, when are real, real separable coordi- 
nates are well-defined, the proposed approach allows to use the same 
notation also when the variables are complex conjugate. 
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The Jacobian matrix and its inverse J arc 



(dz v \ 
\dqi* J 



1 + 



-1 + 



1 - 



V y/l + 2(x+t) 



-1 - 



y^l + 2{x+t) y/1 + 2(x + t) 

y/l + 2{x + t) ) 
1/ 1 + ^ l + 2(x+7) l-y/ l + 2(x + j) 
4 V -1 + V 1 + 2 (^ + *) -1- Vl + 2(x + t) 

hence, the momenta p z and pt are given by 



p x = P 1 +P 2 + 2 P \ P l Pt 

ry± 



Z x — 



Z 1 — Z' 



and the Hamiltonians by 
H = 4 



Pi 2 - P 2 2 



z 1 — z 2 



K = 4 



Z 2 P! 2 - ^P 2 2 



Z 2 — Z 1 



The two Hamilton-Jacobi equations become 



dW 

dW\ 2 
dz 1 J 



i fdW\ 2 _ 1 



dz 2 ) I 



2 z 1 ) 



and by solving them one has (outside the singular line 2(x+t) + l = 0) 
the two equations 



dW 

dz» 



z^h-k 



that admit a separated solution W = Wi(z 1 ) + W2(z 2 ). The two 
functions both satisfy the equation 



zh — k 



that is, by introducing the variable ( = zh — k (assumed not belonging 
to the real negative semi-axis for monodromy reasons), the equation 

dC y ' 2h V ^ 
The solution of this equation is 
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then, if £^ = z^h — k are both real one has 



W = SiC 1 ) + S(C 2 ) = ^(z^h-kf + ^{z*h-k)\ 
while if they are complex conjugate one obtains 

W = l h Re(?)i. 

However, it is not strictly necessary to solve the differential equation 
in order to apply the Jacobi method: 

aw awae owae_ . . 2 2 

aw dwdc 1 dwd( 2 _ x 2 
° k -^k - W~dk W~dk (c } (c) 

where 9 k is constant while 9h = t — tq. 

After substituting the values of £^ as functions of x and t and 
solving with respect to these two variables one obtains the straight 
lines 

6hB k 3 k 2h 4 

1-2^2 h9 2 k + l k ! 



Finally, with the choice h = —k = 2 one obtains the solution 



1 



* = o, x = --(e k e h + i) 

that is a geodesic defined both in the region where z^ are complex- 
valued and in the region where they are real-valued. 



6 Web associated to complex variables 

In this section the level curve of real and imaginary part of complex 
variables associated to Killing tensor in the Minkowski plane are plot- 
ted and compared with the webs corresponding to real eigenvalues. 
The names of coordinate systems are taken from 
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Figure 1: Real and complex webs for SC3 




7 Final remarks 

In this paper a method is presented to obtain, using a set of com- 
plex variables, a complete solution of the Hamilton- J acobi equation 
on a pseudo-Riemannian manifold. The complex variables arise from 
Killing tensors with pointwise distinct eigenvalues and normal eigen- 
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Figure 4: Real and complex webs for SC7 




Figure 5: Real and complex webs for SC8 



vectors, giving a more general definition of characteristic Killing ten- 
sor, starting with a real Hamiltonan is always possible to obtain a real 
complete integral, separated in the complex variables but not in their 
real and imaginary parts. The method can be extended to Hamiltoni- 
ans with scalar and vector potentials and to the study of Schrodinger 
equation through the techniques of j^j and [3]. 
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